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ABSTRACT 

Electron work functions, surface potentials, and electron number density 
distributions and electric fields in the surface region of 26 metals are calculated 
from first principles within the free electron model. The number of "free" 
electrons per atom is taken as the group number as listed in the periodic table. 
Grain orientation effects are not considered. The calculation proceeds from 
an expression of the total energy as a functional of the electron number density 
including exchange and correlation energies as well as a first inhomogeneity 
term. The self-consistent solution is then obtained via a variational procedure 
akin to the Rltz method. Surface barriers are found, in most cases, to be due 
principally to many -body effects, but dipole barriers are small only for a number 
of alkali metals, becoming quite large for the transition metals. As one might 
expect, surface energies are found to be Inadequately described by this model 
which neglects atomistic effects. Considering the simplicity of the model, 
reasonable results are obtained for electron work functions and surface potential 


TM X-52647 



2 


characteristics for all metals studied, maximum electron densities varying by 
a factor of over 60. 

INTRODUCTION 

The wealth of experimental data available today on electronic work functions 
of bare metal surfaces * is not at all matched by theoretical calculations. There 
have been numerous empirical correlations made relating the electron work 
function of metals to atomic volume, compressibility, the first atomic ioniza- 
tion potential, the energy of the lattice, surface energy, and electronegativity. 

2 3 

These efforts are enumerated by Samsonov et al (see also L. N. Dobretsov 

et al and D. Steiner^ et al). Also, some efforts have been made toward formu- 

lating a first principles description of various aspects of this quantity for 

certain metals. However, such calculations of the total (bulk plus surface 

contribution) electron work function have been provided only for the alkali 

metals. The most sophisticated of these is that formulated by Bardeen*® for 

Na. A free electron model was used and the Hartree-Fock equations were 

solved approximately. 

This is in contrast to the progress made in overlapping areas. For example, 
11 12 

many -electron and atomistic effects have been included in theoretical 

studies of bulk metallic properties of many metals. Likewise, many -electron 

effects and some atomistic effects have been included in the theory of adsorption 
13 

on metals using modern formulations of the many -electron problem. 

A second topic considered here which is related to the electron work function 
is that of the surface potential. Recently, a calculation of the surface poten- 
tial of Na which refines Bardeen's work by making use of a modern many -electron 

11 14 

formulation has been provided by Loucks and Cutler (see also Ref. 15). 
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However, these authors neglect the effect of the surface dipole potential and 

place an infinitely high potential barrier at the surface in order to calculate 

wave functions. The first assumption may well be reasonable for Na, but 

it will be shown that dipole barriers cannot be neglected for most of the metals 

studied here. The second assumption, of course, rules out self-consistency. 

1 fi 17 

More recently, Bennett and Duke ’ ‘ have introduced self-consistency into 

a many -electron calculation of the one -electron potential at a bi-metallic interface. 

A small step is made here toward bringing bare surface work function 

theory up to the level of sophistication of neighboring fields, and in the process 

to gain a greater knowledge about metal surface properties in general. A 

calculation of the work function is presented here for 26 metals including Na 

using the jellium model. In addition the electrostatic (double layer) barrier, 

representative electric fields, electron number density distributions and one- 

electron potentials were calculated for the surface region. The jellium or free- 

electron model is used here so that many surface parameters can be calculated 

rather simply. Conclusions can then be made as to which surface characteristics 

are adequately described in this model and which require further sophistication 

in their description. Also our understanding of the metal surface can be considerably 

18 

enhanced without undue effort. A recent formulation of the inhomogeneous 
electron gas which includes coulomb correlations was used in an approximate 
self-consistent first-principles solution of the model. The number of "free" 
electrons per atom was taken as the group number as listed in the periodic 
table. Grain orientation effects were not considered. 

We found that exchange and correlation potentials make up the major part 
of the surface barrier for most of the metals considered. However, the ordi- 
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nary coulomb potential barriers are significant for all of these metals except 
Cs, Rb, K, and Na. Also, the results obtained using the simple model described 
previously show encouraging agreement with available experimental data for all 
the metals considered. 

This paper is divided into four major sections. In Sec. n. the basic equa- 
tions used are derived. Sec. in. is devoted to a comparison of some of the 
results obtained with existing theoretical findings. Results for 26 metals and 
comparison with experimental data are presented in Sec. IV: Concluding 
remarks are given in Sec. V. 

H. DERIVATIONS 

5 17 19-23 

Following Bardeen, the free electron or "jellium" * * model with 

planar surface (see Fig. 1) is used. Bardeen's use of the Hartree-Fock 
equations is not followed, however, because it provided much numerical diffi- 
culty. Also, since the Hartree-Fock equations neglect antiparallel spin corre- 
lations, attempts to take such correlations into account are necessarily ad hoc 
24 18 

in nature . Hohenberg and Kohn (see also Refs. 25-27) have recently derived 
a powerful formulation of the many-electron problem. This scheme, which uses 
the electron number density as the basic variable, provides considerable simpli- 
fication and includes all many-electron effects in the original formulation. Thus 
it will be used here. 

18 

Hohenberg and Kohn (HK) have shown that the ground -state energy, 

E y , of a confined interacting inhomogeneous electron gas can be written as a 
functional of the electron number density n (r9. Further, they have shown that 
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E y fn] assumes a minimum value for the correct n(j£), if admissible density 
functions conserve the total number of electrons. Thus, n(r) can be deter- 
mined from 



- f N ] = 


=. o 


) 


( 2 . 1 ) 


op 

where fi is a Lagrange multiplier such that fi = dEy/dN, 
and N = l n(jr) djc. 

^ n a 


HK write 


29 


(2. 2a) 

where v(r] is a static external potential, g 0Q £ T sC n ] + E xcC n D» T sC n D is 
the kinetic energy of a system of noninteracting electrons with the same density 
n(r), and E xc [nJ is then the exchange and correlation energy of an interacting 
system. 

HK derive an expansion of GQnJ originally for the case of slowly varying n 

in successive orders of the gradient operator acting on n(r) which can be 

30 

written as follows : 


GC-'J f (4)*f'*** c ^ - 
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(2. 2b) 
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The integrands of the first through third terms on the RHS of Eq. (2. 2b) 
represent respectively the kinetic, exchange, and correlation energy densities 
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of a uniform electron gas of density n. The Wigner interpolation formula was 
used to represent the correlation energy of a homogeneous electron gas at 
metallic densities. The fourth term is the first of the inhomogeneity terms, i. e. , 
those terms containing one or higher orders of the gradient operator acting 
on n. 

Several comments about Eq. (2. 2b)are in order. First, it is shown else- 
31 

where that at least some inhomogeneity terms must be included in a work 
function calculation. It is well known that a simple Thomas Fermi theory 

31 

predicts that the work function of any physical system is zero. We have shown 
that including the homogeneous electron gas exchange and correlation energy 
terms, but not including inhomogeneity terms, leads to a predicted work function 
which is nonzero but is the same for essentially any system. It will be seen 
in the following that the addition of the first inhomogeneity term alleviates this 
anomaly. 

Secondly, the random -phase approximation was used by HK to derive the 

factor l/72n in the first inhomogeneity term. Although the RPA has exhibited 

failings at electron densities as low as those found in conduction bands **, this 

inhomogeneity correction to the total energy apparently has a rather wide 

32 33 

range of applicability as shown by the successes of Kirzhnits and Kalitkin . 
Kirzhnits considered isolated noble gas atoms and Kalitkin compared his 
results with experimental bulk properties of solids. Also the RPA has been 

1 Q Q4 1 C 

used with some success in metal surface theory ’ ’ . Thus it is used here. 

Third, HK note that a ’’gradient” expansion of which the sum of the inte- 
grands in Eq. (2. 2b) is an example does not converge for actual electronic 
systems due to number density variations with position. However, they expect 
it to be useful in the sense of asymptotic convergence for sufficiently slowly 
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varying number densities. A formulation based on the ’’gradient” expansion 
has exhibited some successes even for the case of atoms , where the 

■* 

32 

density variation is rather rapid. Also Kirzhnits investigated explicitly the 

convergence of an expansion of E y [n] in successive powers of fi not including 

correlation energies. He calculated E y [nJ for the argon atom, and found 

in an approximate manner "excellent convergence of the approximation process”, 

at least when his first four inhomogeneity terms were included. Finally, HK 

note that quantum density oscillations are not included in the expansion given 

in Eq. (2. 2b). However, it has been reported for a jellium model with planar 

surface that the Friedel oscillations occurring inside the metal are greatly 

diminished by requiring that the surface potential be self -consistent with the 

electron number density distribution. Since a self-consistent calculation is 

done here, they are neglected. Finally, corrections to the Thomas -Fermi 

37 

equation derived by expansion procedures have been shown by Schey et al. 
to be pejorative in many instances. However, they note that expansions of 
the total energy (as we use here), lead to "remarkable improvement.” 

Keeping only the first inhomogeneity term and combining Eqs. (2. 1) and (2. 2), 
one obtains for our model 



2 0 

where, for self consistency, d ip / dZ = 4rr[n + H(-Z)-n2, n + = (positive jellium 
charge density, H(Z) is the Heaviside (step) function, Z is the cartesian 
coordinate taken on an axis normal to the surface, with Z = 0 at the jellium 
surface, tP — 'V~t /£/) ~h C \ > 

r ~ j m -& i J 
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and, in this case, v(r) is the negative of the potential of the ion distribution. 

Note that in the jellium model, n— n + and all derivatives of n— 0 as 
Z— -oo. Also <p e = electron work function*® h -(SEy/aN)^^ = -/l t. Thus one 
obtains ( <p e is described in Fig. 2) : 



where <p(-») represents the value <p asymptotically approaches deep within 

the metal and <p is set equal to zero at large distances from the metal. 

It should be remembered that this is a many- electron calculation. However, 
25 

Kohn and Sham have shown that it is possible, formally, to replace the 
equations of the many -electron problem by an equivalent set of one -electron 
equations. The effective one -electron potential energy is given 
formally by 
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A comparison of Kirzhnits* first inhomogeneity term in his expansion 
in powers of R of the Hartree total energy and that in Eq. (2. 2b) shows that 
they are identical. Thus, the first inhomogeneity term contributes only to 
T g [n] in 016 RPA * So ’ to °(t n ! 2 ) in E v » 


v'°- ^ - ( 0 ^- 


f/L 

. 0 srfc/vv > , Oos<f J 

(.071 + * 


( 2 . 6 ) 
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as given in Eq. (2. 6) is just the potential energy that one would obtain'* 5 * 
for the highest one -electron energy state of a uniform electron gas of density n 
in its ground state. Thus it follows that, at least to 0(|vn|^) in E y , is 
equivalent to the effective potential energy for a state at the top of local Fermi 
distributions. 

can be obtained immediately through Eq. (2. 6) once the many -electron 

problem is solved, and it is exhibited in several of the figures for the reader 

who is interested in a one -electron calculation. 

In order to obtain n, it is certainly simpler to solve Eq. (2. 3) than 

e.g. , a set of Hartree-Fock equations. However, we will simplify the solution 

of Eq. (2. 1) still further. Let us assume that the extremal of Eq. (2. 1) be- 

8 21 38- 

longs, to a good approximation, to the following family of functions 9 9 


/*= - 'TO. ef ) 




* 


2 . ) 

where /3 is a family parameter. 


z 

z 


(2. 7) 


Note that for every value of /3 the family 2. 7 satisfies certain requirements 
of self-consistency. First, n asymptotically approaches n + in the metal interior 
and zero in the vacuum region outside the metal. Secondly, 

n-n + H(-Z)j]dZ = 0. There are no experimental data on n which provide 
a direct test of the validity of the family 2. 7. It will be shown below, however, 
that the results obtained using these simple functions are in at least as good an 
agreement with experiment as could be expected using a flat surfaced jellium model. 


The corresponding coulomb potential is 
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Eq. (2. 1) becomes 
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or equivantly, 
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(2. 9) 


(2. 10) 


and € y £nJ is the energy density, i. e. , E v W-Am*. <7 is the surface energy, 

or the energy necessary to cleave a metal per unit area of new surface 

formed. Thus a is the total energy of the separate pieces after splitting 

minus the total energy of the unsplit block. 

A simple result of analytical manipulations of the terms on the RHS of 

Eq. (2. 10) up to and including the first inhomogeneity term is provided 

below, except for the correlation energy integral over the range -®° s Z ^ 0. 

39 

This last term was easily programmed, and is designated below as 
I(n + )//3. 



11 


This gives 


40 



where the terms in Eq. (2. 11) are given in the same order as those in 
Eq. (2.2), and where a= 2^^^ (. 079)/n + ^^. 

Thus /3 can be determined by combining Eqs. (2. 11) and (2. 9), and this 
result can be used to determine n and </? via Eqs. (2. 7) and (2. 8). With 
these, the quantities <p e and can be determined immediately from 

Eqs. (2. 4) and (2. 6) respectively. 


m. COMPARISON OF RESULTS WITH THOSE 
OF OTHER CALCULATIONS 


Work Function of Na 

Results obtained here, as well as Bardeen's results are listed in 
Table I. Wigner's uncorrected^ interpolation formula was used in this 
instance, so that a more direct comparison could be made with Bardeen's 
work. 

Considering the different approximations made in the two calculations, 
the agreement is quite good. Notice that the work function and coulomb 
barrier are 0. 39 volts higher than Bardeen's results. No decision can be 
made based on the experimental data as to which theoretical value is more 
accurate. This is because first, the value listed is for polycrystalline Na, 
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and secondly, there are inaccuracies even in the knowledge of this value. 

If our results turn out to be more accurate, then the difference may be 

14 

explained by the fact that, as stated by Loucks , only a partially self- 
consistent solution was achieved by Bardeen with respect to the electro- 
static part of the problem since the exchange potentials were chosen at 

the beginning and held fixed throughout. But, as noted earlier, it has 
17 

been reported that the Friedel oscillations inside the metal are greatly 
diminished by self-consistency requirements. Since these oscillations 
lead to a "humping up" of electronic charge inside the metal which lowers 
the dipole moment and work function, their overemphasis could lead to 
values of these quantities which are too low. 

Finally, it is clear that this calculation supports Bardeen's conclusion 
that the surface barrier of Na is due primarily to exchange and polariza- 
tion forces with ordinary electrostatic forces playing a minor role. 

Surface Energy 

41 

The surface energy for Na was calculated by Huntington using Bardeen's 

potential*®. Table I shows that the surface energy of Na calculated here 

agrees rather well with that calculated by Huntington. Neither is in good 

20 2 19 

agreement with the experimental value of 0. 240 joules/m . Herring 
however, has pointed out that it is not "fair" to compare the surface energy 
a of a jellium metal with an actual metal of the same electron density. 

Table II shows values of a for Na, Li, and K. The disagreement with 
experiment is even more pronounced for Li than for Na and, in fact, o 

_3 

goes negative for n + s 13 x 10 Thus further results weren't listed. 

It should be noted that the electron work functions and surface potential 
characteristics depend on the variation of a (e.g. , da/ d/3 in Eq. (2.9)) and 
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not on the value of a itself. Thus, the fact that the surface energy results 

do not agree with experiment does not imply that the results for the work 

functions and surface potentials should not be trusted. 

Na Surface Potential Characteristics 

To compare our for Na (see figs. 3 and 4) with the results of 

14 

Loucks and Cutler (see their fig. 5), one will have to bear in mind that 

our potential pertains to an electron at the top of the Fermi distributions 

whereas they averaged their exchange contribution. Also, as mentioned 

earlier, they neglect the small coulomb contribution. Their potential curves 

are qualitatively similar to V^*\ except that their curves exhibit noticeable 

damped oscillatory behavior in the interior of the metal. As previously 
17 

noted, Bennett et al have concluded that these oscillations are exaggerated 

by lack of self-consistency. Figures 3 and 4 present curves of and 

relative electron number densities, respectively, in the surface region of 

1/3 

the alkali metals as a function of Z/r e , where r_ = (3/4ffn ) ' . 

S 3 T 

IV. RESULTS FOR SELECTED METALS 
Method of Selection 

It seems reasonable that all metals usually regarded** as ”free- 
electron-like” in their bulk properties could be treated within this model. 
Additionally, the surface properties of even the transition metals have been 
described with a certain degree of success within the free electron model. 
Examples of such successful applications are: the Richardson -Dushman 

equation describing thermionic emission, Fowler -Nordheim vacuum field 

42 43 

electron emission theory , plasma oscillation characteristic loss theory , 

44 

and analysis of periodic deviations in the thermionic Schottky effect . Thus, 
those metals which were in some way amenable to analysis using the free- 
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electron model were chosen for consideration and are listed in Tables 
IV and m. 

The characteristics of the metals enter into the model only through the 
quantity n, . Thus values of n_, the number of conduction electrons per 

T a 

atom, must be designated. For all but the simplest metals this choice 
4 45-47 

is not obvious ’ . However, some properties such as Fermi energies 

48 

of many simple metals are well represented on a free -electron model 
using the group number (as listed in the periodic table), for n The group 

a 

number will be used for n for all metals considered here. It will be 

3 . 

seen that this convention yields surface barrier heights which are consistent 
with experiment. 

For purposes of discussion, the metals are grouped according to common 
properties. The alkali metals, the refractory transition metals, and the 
noble metals are obvious groupings. The rest of the metals can easily be 
grouped according to group number. 

As is seen from figures 5 and 6 the metals considered cover a wide 
range of electron densities, thus, providing a stern test of model and method. 

Electron Work Functions 

Table HI compares our results with the experimental values for poly- 
crystalline metals recommended by Fomenko*. It should be noted that 
there is considerable scatter in the data he collected. 

Several comments are in order concerning the findings listed in Table III. 
First, the theoretical values of <p listed increase with increasing n . 
Secondly, the ordering within groups by experimental work function (e. g. , 
low to high), is generally the same as the analagous ordering by theoretical 
work function. Also, the ordering of groups by average experimental and 



15 


theoretical work functions respectively yields identical results, with the 
exception of the noble metals. Finally, it is seen that theoretical work 
functions of the low n^ metals (principally the alkali metals), are higher 
than the experimental work functions. But for the rest of the metals, the 
theoretical value slips below the experimental value, with the difference 
showing some tendency to increase with n + , again with the exception of 
the noble metals. 

One might be tempted to ascribe the exceptions found in the case of 
the noble metals to the choice of n = 1. That is, although this choice might 

a 

be useful for calculation of certain bulk properties, it may be argued that 

49 44 50 

their surface band structure can be significantly different ’ from 

that of the bulk. However, a recent surface experimental determination 

51 

of the inner potential of Cu gives a value which is consistent with the use 
of n = 1. Inclusion of grain orientation effects may clarify matters. 

A decision based on comparison of experimental data and theory should 
be made as to the accuracy of the jellium model in the prediction of electron 
work functions. This decision is complicated by the fact that there are, 
of course, errors in the experimental data and that grain orientation effects 
are not included in the calculation. From the preceding discussion we have 
seen that there is a general agreement in the ordering of the theoretical 
and experimental work functions within groups and in the ordering of group 
average work functions Further, the deviation of the theoretical work 
functions above or below the experimental values listed is within the range 
of variation conceivably caused by grain orientation effects for the bulk of 
the metals considered. But the entire range of experimental work functions 
is only about 2. 5 volts. Thus, although the theoretical values generally 
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pass the test of comparison with experiment, it is not as stringent a test 
as one might desire. 

But the surface potential characteristics can be turned to for further 
testing. It will be seen in the next section that experimental barrier heights 
vary by about 25 volts. This should provide a much more difficult test 
for the theory. 

Surface Potential Characteristics 

The results for electric field, barrier height, and electrical double 
layer are listed in Table IV. Sample plots of are given in Figs. 3 

and 6-8. Included also on some of the plots is the function -1/4Z. Although 
all surface potentials must asymptotically approach the image potential at 
large distances from the metal, an ambiguity arises because it is not clear 
where to place the Z = 0 plane (appropriate to the function -1/4Z), 
with respect to the jellium surface. Thus the function -1/4Z is not neces- 
sarily the image potential, but can be used for scaling purposes. 

Several trends can be inferred from the results. First, the listed barrier 
heights (maximum value of V^) increase with increasing n + . Secondly, 
although generally the better part of the surface barriers are due to many- 
body effects, the ordinary electrostatic contribution to the barrier is small 
only for the alkali metals through Na. In fact, for some of the refractory 
transition metals, the dipole barrier is more than half of the total barrier. 

A comparison of calculated total barrier heights with experiment for 
electrons at the Fermi level provides another check on the validity of using 
the group number for n . Since the barrier height is quite sensitive to n, 
and since it was only desired to check reasonableness in the choice of n ft , 
listing of experimental values was not made exhaustive or necessarily latest- 
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word. In comparing our theoretical values of barrier heights with experiment, 
one must remember that the effective potential seen by an electron depends 
on its velocity. As previously pointed out V^ applies to electrons at the 
Fermi level and thus values of the surface barrier obtained from say electron 
interference microscopy may well not be descriptive of the maximum mag- 
nitude of V^. Also, as mentioned earlier, it is not necessarily true that 
the experimental barrier height should be given by the experimental value 

of the work function added to the experimental or theoretical bulk Fermi 

44 

energy. For example, D’Haenens and Coomes point out that, following 

this procedure, one would obtain lower total barriers than their (surface) 

experimental values indicate (see, however, Ref. 52). These authors explain 

that the energy -level system could understandably undergo modification at 
20 

the surface Thus wherever (surface) experimental values of the surface 
barrier for electrons near the Fermi level were known to differ significantly 
from the sum of the bulk Fermi energy and electron work function, the re- 
sult of the surface experiment was used in Table IV. A comparison of 
theoretical and experimental barrier heights listed in Table IV shows that 
the values generally agree within experimental error. This lends support 
to the use of the group number for n 

21 

Additionally, a comparison of plasma oscillation theory results with 
the data obtained in surface characteristic loss experiments can be used 
to determine n (see, e.g. , Ref. 43 or Ref. 53). The results of these 

a 

authors support the use of the group number for n for many metals. 

a 

Finally, electric fields were calculated. It follows from Eq. (2.8) 
and the values of /j listed in Table II that the electric field (d</?/dZ) varies 
rapidly with position, always pointing out of the metal. Now in a real 
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metal there are very strong fields in the ion cores (experienced generally 
by the core and not the conduction electrons), which are not present in the 
jellium model. Therefore the electric field is calculated at a somewhat 
arbitrary poir.i outside the metal surface (Z=3(a Q )) where the result should 
be free of strong core field effects. 

The listed values of the fields calculated at the aforementioned point 
increase with n^, increasing by roughly a factor of fifty in going from the 
alkali metals to the refractory transition metals. 

Semiempirical calculations of electric fields as seen by adsorbed 

54 55 56 

particles on molybdenum and tungsten ’ agree rather well with the 
theoretical values obtained here. 

V. CONCLUDING REMARKS 

The following generalizations can be inferred from the results obtained 
here: 

(1) There is approximate agreement between the experimental data and 
the work functions and surface potential characteristics obtained here using 
the free electron model. This lends support to the 

premise that it may be possible to calculate rather accurate values for some 
metal surface characteristics via introduction of refinements to this simple 
model. This may even be so for some of those metals whose bulk character- 
istics are not so easily described, e. g. , the refractory transition metals. 

(2) Many -body effects were found to be of importance in all cases and 
ordinary electrostatic effects are quite strong for many of the metals con- 


sidered. 
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FOOTNOTES 

* Based on a dissertation submitted on April 26, 1968, to Ohio State 
University in partial fulfillment of the requirements of the doctorate 
in Physics. 
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TABLE I. - A COMPARISON WITH BARDEEN’S RESULTS 
FOR WORK FUNCTION, COULOMB BARRIER 
AND SURFACE ENERGY FOR Na 



Double layer moment, 

Work function, 

Surface energy 



eV 


eV 


o 

joules/m 


Here 

Bardeen 

Here c 

Bardeen 

Experi- 

mental 3 ' 

Here 

Bardeen* 5 


0. 786 

0.4 

2.74 

2.35 

2.35 

0. 112 

0. 088 

Neglecting** 

correlation 

0. 978 

~1 






energies 

Neglecting exchange 
and correlation 

3. 12 

~4 






energies 









Q 1 

A value listed for work function is Fomenko's recommended value. 

k The actual calculation of the surface energy using Bardeen's results was 

41 

done by Huntington. 

c 10 

Wigner's uncorrected interpolation formula was used here since Bardeen 

used it. 

d 32 

The h expansion of Kirzhnits was used. 
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TABLE H. - SURFACE ENERGIES 
OF K, Li, AND Na - 
COMPARISON OF RESULTS 


Metal 

o 

Surface Energy (joules/m ) 


Theory 

Experimental 3, 

K 

0. 0688 

0. 146 

Na 

0. Ill 

0.240 

Li 

0. 132 

0. 510 


a 


Ref. 20 



TABLE HI. - WORK FUNCTIONS OF 
SELECTED METALS - 
COMPARISON OF RESULTS 


n + (xi0 3 ) 

(a.u.) 


Ca 2 
Mg 2 
Cd 2 
Zn 2 


1.95 1.32 
3.77 1.27 
6.92 1.24 


8.73 

8.80 

12.6 

6.90 

12.8 

13.8 
19.5 

35.8 


In 3 17.0 

Ga 3 22. 3 

A1 3 26.9 


Work Function (eV) 

Theory 

Experimental 2 

2.64 

1.81 

2.71 

2. 16 

2.76 

2.22 

2.93 

2.35 

3. 11 

2.38 


1.23 3.19 

1.23 3.19 

1.23 3.32 

1.24 3.11 2.80 

1.22 3.33 3.64 

1.22 3.36 4.1 

1.22 3.50 4.24 

1.26 3.75 3.92 


3.30 3.3 

3.40 3.7 

3.44 3.8 

3. 56 3. 96 

3.64 4.25 


3.45 
3. 50 


41.3 
41.6 

56.2 

57.4 

70.4 

84.2 


3.80 

3.81 

3.91 

3.92 
3.98 
4.02 


4. 12 
3.99 
4. 5 

4.3 
5.0 

5.3 


Value listed is Fomenko* s recommended value 
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TABLE IV. - SURFACE POTENTIAL CHARACTERISTICS OF 


SELECTED METALS 1 


Metal 

Electric field b (V/m) 

Barrier height (eV) 

Double layer (eV) 



Theory 0 

Experimental 


Cs 

6.03X10 7 

4.22 


0.258 

Rb 

7. 72xl0 7 

4.54 


0. 327 

K 

9.22X10 7 

4.79 

4. 12 d 

0. 386 

Na 

2. 09x10® 

6.04 

4. 85 d 

0.794 

Li 

4.44x10® 

7.84 

6. 08 d 

1.55 

Ag 

5.62x10® 

8.70 


1.99 

Au 

5.88X10® 

8.74 


2.01 

Cu 

8.62X10® 

10.4 

11. 4 d 

2.91 

Ca 

4.42x10® 

7.82 

5.9 d 

1.55 

Mg 

8.80:10® 

10. 5 

9. 84 d 

2.95 

Cd 

9.48x10® 

10.8 


3. 18 

Zn 

1. 31x10® 

12.9 

15. 24 d 

4.45 

Be 

2. 12X10 9 

17.9 

17. 72 d 

7.76 

La 

8.22X10® 

10. 1 


2. 77 

T1 

1. 05xl0 9 

11. 5 


3.55 

In 

1. 16X10® 

12.0 


3.90 

Ga 

1.47xl0 9 

13.9 


5.06 

AI 

1.71x10® 

15.3 

16. 05 d 

6.00 

Sn 

1. 18X10® 

t 

12.2 


3.98 

Pb 

1. 31x10® 

12.9 

13. 78 e 

4.43 

Ta 

2.34x10® 

19.4 


8.78 

Nb 

2.3010® 

19. 5 


8.84 

W 

2.8010® 

23.0 

22-24 1 

11.4 

Mo 

2.89x10® 

23.3 

11.6 

Re 

3. 25x10® 

26.2 

27-29 1 

13.8 

Ir 

3.6010® 

29.0 


15.9 


a The quantities listed here are obtained self -consistently with 
those listed in Table III. 

^Evaluated at Z = 3(a Q ) 
c Maximum magnitude of 

d Obtained by adding Fomenko's 1 recommended work function to X-ray emission 
band widths as listed In Wilson 46 

*The Fermi Energy as given by Anderson 57 , et al Is added to Fomenko’s 1 
recommended work function 
f Ref. 44 


Density 
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Figure 1. * Electronic and positive charge densities lor the 
jelllum model. 



Figure 2 . - Relation between the electron work function, ^ e . 
the Fermi energy, E fi and the effective one-electron potential 
energy for a state at the top of the Fermi distributions, V< n 
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Figure 4. - Relative electron number density distribution in the 
surface region for the alkali metals. 



Figure 5. - Electron number density 
surface region for W, Al, and Cs. 
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Figure 7. - V* 1 * in the surface region of Al and Mg. 





Figure 8. - in the surface region for selected refractory 



